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Abstract 


The  problem  of  the  discrete  time  detection 
of  a known  constant  signal  in  additive  white  sta- 
tionary Laplace  itoise  is  considered.  Both  the 
Neyman-Pearson  optimal  detector  and  the  linear 
detector  are  treated.  Convenient  expressions  for 
the  distribution  functions  of  the  test  statistics 
are  given.  These  expressions  allow  one  to  deter- 
mine the  performance  of  the  detectors  in  a compu- 
tationally convenient  manner. 


Introduction 


Notice  that  Laplace  noise  has  the  "heavy  tail" 
behavior  associated  with  impulsive  noise. 

The  Laplace  distribution  is  popular  in  sta- 
tistics and  many  of  its  properties  have  been 
studied  [3].  Furthermore,  it  is  used  as  a noise 
model  in  engineering  studies.  For  example. 

Miller  and  Thomas  [IJ  used  Laplace  noise  in  a 
numerical  study  of  relative  efficiency.  Bern- 
stein, et.al.  [4]  conment  on  the  non-Gaussian 
nature  of  ELF  atmospheric  noise,  and  they  give  a 
plot  of  a typical  experimentally  determined  pdf 
associated  with  such  noise  [4,  Figure  10].  This 
experimentally  determined  pdf  is  similar  to  a 
Laplace  pdf,  and  on  a linear  graph  the  difference 
is  barely  distinguishable.  Mertz  [5]  proposed 
the  following  pdf  for  the  amplitude  distribution 
of  impulsive  noise: 


f(x)  = vh''(x  + h)' 


Notice  that  if  we  let 


V = b - 1 , 


1 im  f (x)  = ye* , X > 0 
h-x" 


Recently,  there  has  been  considerable  in- 
terest in  the  detection  of  signals  in  non-Gaussian 
noise.  Although  the  assumption  of  Gaussian  noise 
is  frequently  justified,  such  as  in  UHF,  in  other 
cases,  such  as  ELF,  the  assumption  is  dei'initely 
unjustified.  One  form  of  frequently  encountered 
non-Gaussian  noise  is  that  known  as  impulsive 
noise.  Impulsive  noise  is  typically  character- 
ized as  noise  whose  distribution  has  an  associ- 
ated "heavy  tail"  behavior.  That  is,  the  prob- 
ability density  function  (pdf)  approaches  zero 
more  slowly  than  a Gaussian  pdf.  The  references 
in  [1]  and  [2]  give  a summary  of  some  forms  of 
impulsive  noise  and  situations  where  it  arises. 

In  this  paper  we  consider  the  discrete  time  de- 
tection of  a known  constant  signal  in  additive 
white  Laplace  noise.  That  is,  the  pdf  of  the 
noise  is  given  by 


Thus  the  limiting  case  of  the  Mertz  model  for 
the  amplitude  distribution  of  impulsive  noise  is 
identical  to  the  distribution  of  the  amplitude  of 
Laplace  noise.  Kanefsky  and  Thomas  [6]  consid- 
ered a class  of  generalized  Gaussian  noises,  ob- 
tained by  generalizing  the  Gaussian  density  to 
obtain  a variable  rate  of  exponential  decay.  The 
Laplace  distribution  is  within  this  class  of  gen- 
eralized Gaussian  distributions.  Also,  Dutt- 
weiler  and  Messerschmitt  [7]  refer  to  the  Laplace 
distribution  as  a model  for  the  distribution  of 
speech. 

In  the  next  section  we  present  a brief 
summary  of  the  problem.  Then  in  the  following 
section  we  develop  convenient  expressions  for  the 
distribution  functions  of  the  test  statistic  un- 
der both  hypotheses.  Finally,  some  comments  on 
the  linear  detector  are  given  and  several  exam- 
ples are  given  to  illustrate  the  usefulness  of 
the  results. 


Preliminaries 


We  consider  testing  for  the  presence  or 
absence  of  a positive,  constant  signal  s,  in 
additive  Laplace  noise.  We  assume  that  the  noise 
samples  are  statistically  independent.  (A  re- 
stricted receiver  bandwidth  might  cause  this  as- 
sumption to  be  violated.)  The  problem  is  modeled 
as  the  following  hypothesis  testing  problem: 


Ho:  X.  = n. 


H, : X , = s + n . . 


i = 1,2 N 

s > 0 


Based  on  the  observations  ix.j,  i = 1,2,...,N), 


we  are  to  decide  to  announce  that  the  signal  is 
absent  or  present.  The  quantity  a will  denote 
the  probability  of  false  alarm;  that  is,  a is  the 
probability  of  incorrectly  announcing  H.j . Sim- 
ilarly, s,  the  detection  probability,  is  the 
probability  of  correctly  announcing  Hi. 

The  Neyman-Pearson  optimal  detector  is  a 
detector  which,  for  a fixed  o,  will  maximize  B. 
The  optimal  detector  for  our  problem  is  well 
known  [8],  and  is  illustrated  in  Figure  1.  The 
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observations  are  passed  through  a zero  memory 
nonlinearity  g(-),  and  then  suimied.  The  result 
is  then  compared  to  a threshold.  The  threshold 
T is  chosen  to  give  the  desired  false  alarm 
probability.  The  nonlinearity  g(-).  illustrated 
in  Figure  2,  is  given  by  the  following  expres- 
sion: 

ys.  X > s 

g(x)  = 2yx-ys,  0 1 X ^ s (2) 

-ys,  X < 0 

For  the  optimal  detector,  the  test  statis- 
tic t is  given  by  the  following  sum  of  indepen- 
dent, identically  distributed  random  variables; 

N 

t = I g(x.)  . 
i=l  ’ 

If  the  distribution  of  this  sum  were  known,  then 
the  detection  and  false  alarm  probabilities  could 
be  found,  and  the  performance  of  the  detector 
would  be  known.  However,  past  attempts  at  ob- 
taining a simple  expression  for  this  distribution 
have  not  been  successful.  A lengthy  and  complex 
recursion  scheme  for  obtaining  this  distribution 
has  been  considered  by  Miller  and  Thomas  [1,  9]. 
If  N were  sufficiently  large,  the  Central  Limit 
Theorem  would  apply,  and  the  distribution  of  t 
wuld  be  approximately  normal.  However,  the 
small  sample  performance  of  the  detector  would 
still  be  unknown  (see,  for  example,  [1,  10]). 
Alternatively,  one  could  establish  bounds  on  the 
detection  and  false  alarm  probabilities,  and  thus 
establish  a bound  on  detector  performance;  or 
Monte  Carlo  simulation  may  be  employed.  In  gen- 
eral, however,  it  would  be  desirable  to  have  a 
convenient  expression  for  the  probability  dis- 
tribution of  the  test  statistic  t. 

The  Neyman-Pearson  Optimal  Detector 

In  this  section  we  derive  an  expression  for 
the  distribution  of  the  test  statistic  for  the 
Neyman-Pearson  optimal  detector.  The  test  sta- 
tistic is  obtained  by  passing  each  of  the  obser- 
vations through  the  nonlinearity  g(-).  given  by 
Eq.  (2),  and  summing  the  outputs. 

We  first  consider  the  case  of  no  signal, 
i.e.  Hq.  If  Xi  has  a Laplace  distribution  given 
by  Eq.  (1),  then  g(Xi)  will  have  the  following 
distribution  function: 

F(x)  = i u(x+ys)  + if  exp[-  I(ys+v)]G{^  ]dv 
2 4J.„  2 tys 


+ ^ e-''"  u(x-ys)  , (3) 

where  u(-)  denotes  the  unit  step  function  given 
by 

1 , X > 0 

u(x)  = 

0,  X < 0 

and  G(-)  denotes  the  gate  function  given  by 


The  distribution  function  of  the  test 

statistic  t is  given  by 

where  Fi(x)  = F(x). 

The  Fourier-Stieltjes  transform  of  Eq.  (3) 
is  given  by 


F(u.)  = f e'J'*’V(x)  , 

where  j denotes  the  imaginary  unit.  A straight- 
forward calculation  yields  that 

F(u)  = exp|-  ^]^'=oshj^|2’^jujysj 

y . (4) 

1 + 2ju  / 

Letting  Ffj(ui)  denote  the  Fourier-Stieltjes  trans- 
form of  Ffj(x),  we  get  that 

F^tco)  = [no;)]*^  . (5) 

Using  the  binomial  expansion,  it  follows  from 
Eqs.  (4)  and  (5)  that 


Expressing  the  hyperbolic  sine  and  hyperbolic 
cosine  in  terms  of  complex  exponentials  results 
in 


Using  the  binomial  expansion  again  yields 
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-k  „k-N 


f,(.i . I 


^0  IpJ  ^ 

(^'’“’H'’] e»p[[|^j'^] (N-k-q )j  . 


p=( 

N-k 


A straightforward  simplification  results  in  the 
following  expression: 


Ffj(uj)  = I 


k=l 


^lUJ 


,-(N+k) 


I ' (-1)'^ 


piolP 


i>-^exp[-(p>q).s]P-rtj-<'^-^P-^q)rsI 


N-k 

q 


+ 2'^  I ([!^exp{-ii^s)exp[ jiij(N-2m)Ys] 
m=0\"'/ 

Let  A(w)  = F^(ai)  - 2'^ 

• exp[  ju)(N-2m)Ys]  . 

Notice  that  A(o))  belongs  to  Lg  and  thus  possesses 
an  inverse  Fourier  transform  fl(x),  defined  as  a 
limit  in  the  mean. 

• I (^‘'^)exp[-{p+q)Ys] 

q=0  \ 

. _L  [ “ exp(M(S-2p-2q)rS+xV 

Using  contour  integration,  the  residue  theorem, 
and  Jordan's  lemna  [11],  we  obtain,  for  k^l , 


1 

^ J. 


exp{jair(N-2p-2q)YS^x]}  ^ 

[M‘ 


,k-i 


£^x-»(N-2p-2qjYS] exp{-  ^x+(N-2p-2q)Ys]] 

(k-1)!  ^ 

• u[x+(N-2p-2q)YS] 

After  simplification,  we  obtain 


, N ,-{N+k) 

A(x)  = N!  exp[-i(x+N>s)]  [ 

^ k=l  (k-1)! 

. I (-I)**  [x^(N-2p-2q)Ys]‘^-^ 

p=0  p!(k-p)!  q=0  q!(N-k-q)l 


Let 


• u[x+{N-2p-2q)rs] 


B(x)  = f A(v)  dv 

J _LO 

Upon  performing  the  integration 
obtain 


(see  [12]), 


we 


uu  ba  I II 

B(x)  = N!  exp[-l«YS]  I 2'^'^*'''  J -- 

k=l  p=0  p!(k-p)! 

u^x■^(N-2p-2q)YS] 
q=0  ql(N-k-q)! 

_Y  [(N-2p-2q)Ys]'^~'-^ 

r=0  (k-r-1)! 

• I {(-l)'‘-“[(N-2p-2q)Ysf'“ 

w=0  (r-w)!  ' 

• exp[|^(N-2p-2q)rs]  - x^~'^  } . 

Thus  we  get 

AM  = I A(x)e‘J‘'”‘  dx  = j e"J“”‘  dB(x)  • 
Letting 

C(x)  = 2‘^  I exp(-mYs)  u[x-(N-2m)Ys]  , 

m=0 

we  get 


dC(x)  = I exp(-mYs) 

m=0  ' / 

exp[ju)(N-2m)Ys] 


foo  . ftl  N 

e-J“>x  -•  ■ -N  - 


Thus 


F^(u,)  - “e’j'"'' 


anil 
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Therefore,  we  see  that  Ff^(x)  = B(x)  + C(x). 

Letting  F^°^(x)  denote  the  distribution 
function  of  the”test  statistic  under  the  hypo- 
thesis H , we  have,  after  a straightforward  ma- 
nipulation. 


•j^exp[-(p+q)vs]-exp[-^x+N>s)] 

• ®k.i{^'^(N-2p-2q)Y^^J  •u[x+(N-2p-2q)ys] 


(6) 


N N 

+ 2"'*  [ 
m=0 


exp(-n»Ys)-  u[x+(N-2m)Ys]  , 


where  e.(-),  the  incomplete  exponential,  is 
defined^as 


ek(x) 


k m 
m^O^ 


We  now  consider  the  signal  present  case, 
i.e.  . We  let  F|^^^(x)  denote  the  distribution 

function  of  the  test  statistic  under  Hi-  Since 
the  Laplace  pdf  is  symmetric,  it  can  be  shown  [1] 
that 

F^’’(x)  = 1 - F|f°>(-x)  . (7) 

Eqs.  (6)  and  (7)  thus  completely  determine  the 
performance  of  the  Neyman-Pearson  optimal  detec- 
tor. 

The  Linear  Detector 

By  a linear  detector,  we  mean  a scheme  such 
as  that  illustrated  in  Figure  1,  but  where  the 
function  g(-)  is  g(x)  = x.  That  i^,  the  test 
statistic  is  simply  the  sum  of  the  observations. 
The  linear  detector  is  Neynan-Pearion  optimal  for 
Gaussian  noise  and  is  a commonly  used  detector. 

Consider  the  signal  absent  case,  i.e.  H^. 

In  this  situation,  the  test  statistic  is  given 
by 

N 

t = , 

i*l  ’ 

where  the  are  independent  identically  dis- 
tributed random  variables  with  the  pdf  of  Eq.(l). 
Let  Pn(x)  denote  the  pdf  of  t.  Then  we  have 
[3,  p.24] 


G^°’(x)  = 


•] 


• Vk-i('^")l 
1 - Gl°h-x}  , 


X ^ 0 
X < 0 


(8) 


In  the  signal  present  case,  the  test  sta- 
tistic is  given  by 


N 

t = [ X,  + Ns  , 
i=l  ^ 


where,  once  again,  the  Xj  are  independent  iden- 
tically distributed  random  variables  with  the 

density  function  of  Eq.  (1).  Let  denote 

the  distribution  function  of  the  test  statistic 
of  the  linear  detector  under  . Then  we  have 

' gJ||^>(x)  = G^°'(x-Ns)  . (9) 

Eqs.  (8)  and  (9)  completely  determine  the  perfor- 
mance of  the  linear  detector. 

Examples 

As  an  example,  assume  that  y=1,  s=l , N=6, 
and  that  we  set  a=0.I.  Then  for  the  optimal 
detector,  we  find  that  the  threshold  is  given  by 
T=0.4044  and  the  detection  probability  is  given 
by  6=0.809.  For  the  linear  detector,  we  get  that 
the  threshold  is  T-4.3265  and  the  detection  pro- 
bability is  given  by  6=0.696. 

Now  we  set  cf=0.01,  with  y=1  , s=l , N=6,  as 
above.  For  the  optimal  detector,  we  find  that 
T=2.7321  and  that  6=0.418.  For  the  linear  de- 
tector, we  get  that  T=8.43  and  6=0.23. 

As  another  example,  assume  that  y'1,  s=1, 
N=36,  and  that  we  set  a*0.1.  For  the  optimal 
detector,  we  find  that  the  threshold  is  given  by 
T=-6.9234  and  the  detection  probability  is  given 
by  6=0.999957.  For  the  linear  detector,  we  get 
that  T=10.824  and  6=0.9982. 

Now  we  seta=0.01,  with  y=1  , s=l,  and  N=36. 
For  the  optimal  detector,  we  get  T=-1.5093  and 
6=0.998.  For  the  linear  detector,  we  find  that 
T=19.9  and  6=0.97. 

As  above,  we  let  y'1.  s=l , N=36;  and  now 
we  set  a=0.001.  For  the  optimal  detector  we  find 
that  T»2.503  and  6=0.9837.  For  the  linear  de- 
tector, we  get  that  T=26.7935  and  6=0.8625. 
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